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Abstract 

The triple gluon field strength operator represents the only genuinely gluonic CP 
conserving term which can appear at dimension-6 within an effective strong interaction 
Lagrangian. Previous studies of this operator have revealed that its effect on gluon scat- 
tering is surprisingly difficult to detect. In this article, we analyze the impact of upon 
top quark pair production. We find that it will generate observable cross section deviations 
from QCD at the LHC for even relatively small values of its coefficient. Furthermore, G^ 
affects the transverse momentum distribution of the produced top quarks more strongly 
at high energies than dimension-6 four-quark and chromomagnetic moment terms in the 
effective Lagrangian. Top-antitop production at the LHC will therefore provide a sensitive 
and clean probe for the elusive triple gluon field strength operator. 
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1. Introduction 



The quest to uncover signs of new ptiysics beyond tlie Standard Model is being pursued 
across a broad front. One important area in which the search has been underway for many 
years is the strong interaction sector. Quantum chromodynamics has so far passed every 
experimental test and explains all strong interaction phenomenology at currently accessible 
energies. However, new color interactions beyond those of conventional QCD could emerge 
at higher energy scales from a number of different sources. For example, a wide range of 
theories beyond the Standard Model contain novel particles which couple to quarks and 
gluons. Such particles could be squarks and gluinos in supersymmetric theories [|l| , colored 
technihadrons in non-minimal technicolor models 0, or fermions that transform according 
to representations of color SU{3) other than the fundamental triplet in certain electroweak 
symmetry breaking scenarios 0. New strong sector physics might alternatively originate 
from quark or gluon substructure Preon exchange between composite states would 
generally induce nonstandard quark or gluon couplings. In short, the variety of novel color 
interactions which have been proposed in the literature provides many possible sources of 
interesting physics. 

A complete description of any new strong interaction phenomena requires a funda- 
mental theory beyond the Standard Model. However at energies below its characteristic 
scale A, the new physics may be studied within an effective field theory framework. Its low 
energy effects can be reproduced by supplementing the renormalizable terms in the QCD 
Lagrangian with higher dimension operators. The resulting effective Lagrangian 

^eff = C^co + + ^J2cl'^(^')0l'\^)+0{^) (1.1) 

i i 

generally contains all terms consistent with SU{3) x SU{2) x U{1) gauge invariance and 
any imposed discrete or global symmetries. The operators O^l^^ and their dimensionless 
coefficients C^^^ are grouped together in eqn. ( |1 . Ij ) according to their mass dimension n 
and multiplied by a factor of l/A"^"^ so that their overall combined dimension equals 4. 
As a result, the most important nonrenormalizable terms in Ces are those of dimension-6 
since they are least suppressed by inverse powers of the high energy scale A. 

In the quark sector, one can form a large number of dimension-6 operators with 
different chiral and color structures The qualitative impact of these terms upon 
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quark scattering has conventionaUy been assessed by including just one representative 
four-quark operator such as 

OfXu.r^ = l{q,rq.){q,l,q.) (1-2) 

into the effective Lagrangian and setting the magnitude of its coefficient equal to 4n. The 
scale A associated with the four-quark operator has then been constrained by comparing its 
effect upon the theoretical prediction for the inclusive jet cross section with experimental 
measurement. A recent such comparison yields the bound A > 1.4 TeV 0. 

The number of nonrenormalizable terms which arise at dimension-6 in the gluon sector 
is much more limited. One can build only two gauge invariant operators out of covariant 
derivatives L"^ = d^^-igsG^Ta and gluon field strengths G^" = d^'G^^-d^'G^+gJabcG'^G''^ 
that preserve C, P and T [g] : 

Op"* = 9sfabcG'^yGl)^G^^ (1.3a) 
of) = \d^GI,DxG^,\ (1.36) 

Other candidates such as ggdabcG'^^G'j^-^G'^^ or ^D^C^j^DxC^" either vanish or reduce to 
combinations of 0\ and O2 . The triple gluon field strength term in ( |1.3a| ), which we 



shall name G^ for short, represents a true gluonic operator. It contributes for instance 
to the distinctive three-point vertex in the nonabelian color theory. On the other hand, 
the double gluon field strength operator in (|1.36| ), which we will call [DG)'^, is not really 



gluonic. The classical equation of motion 

i^^or = -9s Yl (1.4) 

flavors 

relates its S-matrix elements to those of a color octet four-quark operator H]: 

of ^-^ f i^^^T^^) (^^'^«^)- (1-5) 

flavors 

The double field strength operator thus affects parton processes involving external quarks 
rather than external gluons. While the origin of O2 could significantly differ from those 
of other four-quark operators, its impact upon parton scattering is not very different from 
theirs. So at dimension-6, only one genuinely gluonic term which conserves CP may appear 
within the effective Lagrangian. 
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The list of CP even gluon operators grows at dimension-8 but remains manageable 
in size. It is useful to classify them according to the number of field strengths that they 



contain [10]. Only one independent operator can be built out of two field strengths and 



four covariant derivatives: 



of) = ^D^Gl^D^Dy^G^:". (1.6a) 

Any other candidate in this category such as ^D'^G'^j^D'^G'^^ may be reduced via the 
Jacobi identity DxG'^^^ + D^G"^^ + Dj^G^^^ = to plus operators with more than two 
field strengths. There are two possibilities for dimension-8 operators with three gluon field 
strengths and two covariant derivatives: 



^2^'' - gsfabcG^^DxGl'' D'^G'l^ 
o'i^ = gsfabcG^^GlxD'^Gcii- 



;i.66) 



(1.6c) 



Finally, six independent operators containing four field strengths can be formed [ pT| : 

2 

0(8) _ Qs 
2 

0(8) _ 9s 
5 ~ l^^iiv^a ^\a^h 

2 

0(8) _ 9s 
6 ~ ~^^ii.v^b ^Xa^h 

2 

0(8) _ ^7. 

2 

0(8) _9sr r 
8 — -^O-abeacde^^jy'^b '-^ Xa^ d 

2 

9 — ^"afae"cde<J"^!y<J"b ^ Xa^ d ■ 

The equation of motion in ( |1.4| ) relates the S-matrix elements of 0\ and O2 to those of 
operators with four quark fields. So among all the dimension-8 operators, only 0\ - 0\ 
are genuinely gluonic and affect scattering processes with external gluons. 

Present day experiments at the Fermilab Tevatron are more sensitive to the presence 
of nonrenormalizable quark operators in the strong interaction effective Lagrangian than 



^ Operator O3 does not reduce to plus a four field strength operator as previously 



reported in ref. |llO|] . 
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to their gluonic counterparts. But at higher energy machines such as the LHC, the gluon 
content of coUiding hadrons at smaU parton momenta fractions dominates over that of 
valence and sea quarks. Future coUider experiments wiU therefore probe possible deviations 
from the Standard Model much more sensitively in the gluon rather than quark sectors. In 
this article, we will investigate the prospects for detecting such deviations. In particular, 
we will focus upon searching for the triple gluon field strength operator G^. 

The most obvious channel in which to look for might seem to be gg —>■ gg scatter- 
ing. However, the helicity structure of the amplitude for this process which involves the 
gluon operator is orthogonal to that of pure QCD. They consequently do not interfere at 
0(1/A^) 1^. At 0(1/A^), the dimension-8 operators O^^ - O^^ enter into the gg gg 
amplitude and contaminate any signal from O^^ |T^. So looking for G^ in lowest order 



gluon scattering is difficult. As we shall later see, the triple field strength operator does 
not mix at one-loop order with any four-quark operators. Its presence in the effective 
Lagrangian thus cannot be indirectly inferred from pure quark processes either. 

These obstacles to observing G^ have motivated consideration of other possible detec- 
tion methods. For instance, it may be probed in Z — > 4 jets decay at LEP |[1^,|13[. This 
option is limited, though, by LEP's relatively low energy. Alternatively, one can look for 
the effect of G^ upon gg ggg scattering in three-jet events . The task of identifying 
its signature in this next-to-leading order channel represents a formidable experimental 
challenge. The triple gluon field strength operator has thus remained stubbornly difficult 
to detect. 

In this paper, we take a new approach to searching for the G^ operator which pos- 
sesses significant advantages over those previously explored. Specifically, we investigate 
the effect of the triple field strength term upon top quark pair production. As we shall 
demonstrate, the impact of G^ at LHC energies is sizable for even relatively small values of 
its effective Lagrangian coefficient. It dominates in this mode over four-quark competitors 
with comparable coefficients. Moreover, contamination from gluon terms of dimension-8 
and higher is small. So ti production provides a sensitive and clean probe for the elusive 
G^ operator. 

Our article is organized as follows. In section 2, we discuss the renormalization group 
running of all the dimension-6 operator coefficients that infiuence top quark pair produc- 
tion. We then compute in section 3 the ti differential transverse momentum cross sections 
at the LHC and Tevatron for representative values of their operator coefficients and com- 
pare with corresponding QCD results. In section 4, we quantify the difference between the 
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transverse momentum cross sections calculated with and without the gluonic operators in 
the effective Lagrangian as a function of their coefficients. In section 5, we investigate the 
operators' impact upon ti angular distributions. Finally, we summarize our findings and 
present our conclusions in section 6. 



2. Dimension-6 operator coefficients 

Many forms of physics beyond the Standard Model could give rise to nonrenormaliz- 
able gluonic operators within the effective strong interaction Lagrangian. In the absence of 
a fundamental theory, the coefficients of such operators are a priori unknown. Even when 
a particular high energy model is specified, their values may prove intractably difficult 
to calculate. So in this work, we will simply study the feasibility of detecting the lowest 
dimension operators as a function of their coefficient values. 

Since the dimensionless coefficients and dimensionful scale appear together within the 
effective Lagrangian ( p.. 11 ), they cannot be distinguished at tree level. We are therefore 
free to fix the value for the high energy scale. Motivated by the expectation that a new 
fundamental layer of physics awaits discovery in the TeV regime, we will take it to be 
A = 2 TeV. This relatively large value for the effective theory cutoff ensures that our low 
energy analysis will be valid over almost the entire practical energy range of present and 
anticipated hadron colliders. 

Although we will generally treat the coefficients of the gluon operators as free parame- 
ters, we note that there is an important special case in which they can readily be computed. 
This occurs when the operators are generated through diagrams involving heavy colored 
particle loops like those illustrated in fig. |l|. At energies below the particle's mass, the be- 
havior of such graphs may be reexpressed in terms of local but nonrenormalizable gluonic 
operators. The values for their coefficients are then determined by performing a matching 
computation. For example, if the particles running around the loop in fig. |^ are fermions 
of mass A and belong to representation R of color SU{3), the sum of the diagrams contains 
terms which match onto the dimension-6 operators in eqn. ( p..3| ) with coefficients 

C(6)(A) = 

' 360 TT .2.1) 



^ 15 TT 
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where K{R) = X]a=i Tr[^a(-R)T'a(-R)]/8 denotes the index of R. The corresponding result 
for colored scalars is given by 

^ 720 TT /2 

2 ^ ^ 120 TT ■ 

Unfortunately, the numerical values of these coefficients are too small in any reasonable 
theory to produce detectable deviations from the Standard Model. We will therefore have 
to be content with constraining mechanisms other than perturbative radiative corrections 
which can generate gluonic operators with larger coefficients. 

Whatever values C^^ and C'^^ may assume at the scale A, they must be evolved down 
to the energies at which and are probed using the renormalization group. We 
consequently need to enlarge our operator basis so that it closes under renormalization. It 
is sensible to apply the equations of motion to reduce the operator set as much as possible. 
We then find that 0\ and O2 do not mix with each other under the action of QCD at 

(6) 

one-loop order. Instead, 0\ runs into itself and the chromomagnetic moment operator 



oS') = 9sm,qa^^T<^qGl,. (2.3) 

flavors 



Their evolution is governed by the 2x2 anomalous dimension matrix 

oS"' of' 

..-"^.i'"' _ ° .\A^Oi,t, (2.4) 



Of> \ 9/2 7 + 2n//3 / S^r^ 



where n/ denotes the number of active quark flavors [lIl|JT5|] . In the sector, the 

equations of motion convert the double gluon fleld strength operator into the color octet 
four-quark term in ( |1 . 5| ) . then mixes at one-loop order with its four-quark counterparts 



flavors 

of = f E (^7m9)(97^9) (2.5) 



flavors 

,2 



of = I E (^7.7^?) {qi'l'q) 

flavors 
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through the 4x4 matrix 





o'i" 






Of' 




^(6) 


/ 311/36 -2n//3 


5/4 





2/3 \ 






41/36 


35/4-2n//3 


2/3 





9l 


of) 


4/3 


6 


11 -2n//3 





87r2 


of 


V 22/3 








11 - 2n//3 / 





7/1= Jfi^ ' „ „ fz5+0{gl). 

(2.6) 

The coefficients of operators - satisfy the integrated renormahzation group 



equation 



E|- 7^(^s 
exp / dgs ^ 



I3{9s 



(2.7) 



where /3{gs) denotes the QCD beta function, i The leading log running of these coefficients 
can be substantial. For example, if they assume the values 



(2.^ 



{Cf^ , ) , ) , ) , ) , ) ) ( A) = (1 , 1 , 1 , 0, 0, 0) 
at the canonical A = 2 TeV scale, the coefficients run down to 

(0.7858, 0.7458, 0.8856, -0.0294, 0.0003, -0.0152) 

at the top-antitop threshold, i Renormahzation group evolution therefore generally sup- 
presses the operator coefficients. 



(2.9) 



3. Top quark pair production 

The sensitivity of a hadron collider experiment to nonrenormalizable terms in the 
effective strong interaction Lagrangian is a function of the accelerator's center-of-mass 

^ We take the constant of integration which enters into the integrated QCD beta function to be 
as{Mz) = 0.125 [ p!7[ . The corresponding value for the strong interaction fine structure constant 
at A = 2 TeV is as (A) = 0.086 which implies ^^(A) = 1.04. Since this coupling is very close 
to unity, the value of C[^^ (A) is essentially independent of the number of powers of gs that are 
included into the operator's definition. 

^ We adopt the recently reported CDF central value mt = 174 GeV for the mass of the top 
quark fl^ . 
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energy ^/s and luminosity C. The Fermilab Tevatron currently operates at ^/s = 1.8 TeV 
and has collected approximately 20 pb~^ of data. These values for and / Cdt are 
too low to allow for any significant limit to be placed upon the coefficient of the triple 
gluon field strength operator G^. However, the prospects for probing will substantially 
improve with the advent of the LHC. This machine is projected to run at ^/s = 14 TeV 
and collect approximately 30 fb~^ of data per year. At such high energies, the gluon 
content of the colliding hadrons at small parton momenta fractions dominates over that 
of all other partons. Sensitivity to the triple field strength operator will consequently be 
greatly enhanced. 

As mentioned in the introduction, looking for G^ in gg gg scattering is problematic. 
The only other 2 — > 2 parton process that maximally benefits from the large numbers of 
gluons in a high ^/s initial state is gg qq. The 0(1/A^) interference between the QCD 
and G^ amplitudes for this mode is proportional to the squared quark mass m^^. It is 
therefore greatest when the final state quark fiavor is top. The top quark is also the 
easiest to tag since its leptonic decay channel can produce a high energy, isolated lepton 
in conjunction with a bottom quark. This distinctive signature cuts down on genuine 
backgrounds as well as false identifications. If in the future bb and even cc final states can 
be positively identified, then the signal for G^ will only be enhanced and all our results 
can be easily applied to their study. But for now, we shall examine the impact of the 
nonrenormalizable terms in Ces upon just top quark pair production. As we shall see, this 
mode provides a sensitive and clean probe for the G^ operator. 

Top quark pair production proceeds at tree level through the parton reactions gg tt 
and qq ti. We first consider the gluon channel. The lowest order QCD graphs that 
mediate gg ti scattering are illustrated in fig. At 0(1/A^), the chromomagnetic 
moment and triple gluon field strength operators oj^^ and 0[^^ contribute through the 
diagrams shown in fig. |]b and fig. ^jc. They also enter into the gg —>■ ti amplitude at 
0(1/A*) via the double operator insertion graphs displayed in fig. |^. The only other 
possible dimension-6 operator tree diagrams contain insertions of the double field strength 
operator O2 ■ But the sum of such (DG)'^ graphs vanishes as guaranteed by the S-matrix 
relation in ( [1.51 ). 
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After adding togetlier all the diagrams in fig. ^ and squaring the total gg — > tt ampli- 
tude, we find @ 



J] \A{gg^tt)\' 



3 (m 
4 



ml 



u 



) 1 m1{s — Ami) 



1 

+ 6L 



24 {ml - t){ml 



u 



iu — m1(3i + u) 



mt 



iu 



2m\i mt'^ 



+ 



s{m\ - 1) 



+ 



iu 



tu — mfit + 3w) — mt 
{mf - 

2mfu + mt^ 



u 



u 



1 



1 ^(6) m| {4s^ - 9iu - 9mp + 9mt'^) 
3 ° (mf — i) (mf — u) 

9 



(6) 



+ 



A4 



m1{i - -u)^ 
(m^ — i) {ml — u) 
^Q)2m'i{Usiu + m2(31s2 - 36M) - SOmt^s + 36mt6) 



{m1-t){m1-u) 



9 -^S 

m^s 



U) 



(3.1) 

The bar appearing over the summation symbol on the LHS of ( |3.1| ) implies that the squared 
matrix element is averaged (summed) over initial (final) spins and colors, while the prime 
indicates that \A{gg tt)p is divided by g\. Notice that all of the nonrenormalizable 
operator terms except the last one are proportional to m^. The corresponding terms in 
the squared gg — gg amplitude for quark fiavors other than top are therefore negligible 
by comparison. The last term in ( |3.1j ) is significantly enhanced by its 27/4 prefactor. 
Moreover, it increases quadratically in the partonic Mandelstam invariants s, t and u. In 
contrast, the other 0(1/A^) terms only grow linearly, while the (9(1/ A^) terms approach 

- (6) 

a constant. So away from the tt threshold and over large regions of C\ parameter space, 
the oj^"* operator's squared amplitude is much larger than its interference with QCD. As 
a result, the impact of upon ti production depends mainly upon the magnitude rather 
than sign of its coefficient. 

( a\ 2 _____ 

Since the gluonic term proportional to C\ in ( |3.1j ) is surprisingly large, one may 
question whether other 0(1/A'^) terms arising from dimension-8 gluon operators could 



The results displayed here in eqns. (3.1) and (3.3) correct errors in eqns. (14a) and (14b) of 



ref. [18| and eqns. (4.4) and (4.5) of ref. ||10[]. The pure QCD terms in these formulae agree with 
the results of ref. lit 
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be significant as well. The answer is generally no. Recall that among the dimension-8 
operators listed in eqn. ( p..6|) , O^^ and 02^'' enter at tree level into processes involving at 
least four quarks, while vertices from 0\ - Og contain at least four gluons. So only O3 
can affect gg tt scattering at lowest order. The interference between its amplitude and 
that of pure QCD yields the 0{1/K.'^) term 

gMto^g)|'^...-|^ "f'f-f.. . (3.2) 
^-^ 8 [mf — t)[mf — u) 

This dimension-8 term has a much smaller prefactor and increases more slowly with s, t 
and u than its dimension-6 competitor. So unless its coefficient C3 is more than an order 
of magnitude larger than the operator is not likely to obscure any signal from 

of\ 

We now turn to consider the quark process qq — ti. The chromomagnetic moment, 
double gluon field strength and four-quark operators in our basis contribute at lowest order 
in the strong interaction coupling to this channel as indicated in fig. |^. Neglecting all quark 
masses except that of the top in the amplitude sum, we find 

\Am ^ tt)? = ^ [^'^ + + 4m?s - 2mt^] 

+ ^ [ACl^^mls + ) {P + ix' + Amis - 2m/) + cfs{i - u 



4 

+ 



9A4 



9sA2 

2 

8Ci^) mlitu + 2mls - mt'^) / s 
+ ^Cl^^cf'^mls + %cl^^cfml{i - u) 

+ (CfVicf')(F + ^.^-2m,^) 



+ (2Cf)C;i')+cf)cf))s(t 



u 



(3.3) 

Unlike the gluonic scattering result in (|3.1|) , this expression contains no anomalously large 



0(1/A^) term. So we expect that the effect of dimension-8 and higher operators upon 
qq tt scattering is small. 

The squared amplitudes in ( p.l|) and ( p.3|) enter into the partonic differential cross 
section 

^-^^^!^ = ^^\A{ab^tt)f (3.4) 
dt 
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which appears in the fuU hadronic differential cross section for top-antitop production |2C]: 



d^a , -^ ■r-^ , ^ r daiah ^ tt) 

[AB tt) = 2p± 2^XaXbfa/A{Xa)Jb/B{Xb) ■ (3.5) 



dy3dy4dp± ^ dt 

ab 

The partonic cross section is folded together with distribution functions fa/Ai^a) and 
fb/sixb) that specify the probability of finding partons a and b inside hadrons A and B 
carrying momentum fractions Xa and Xb- The product is then summed over initial parton 
configurations. The resulting hadronic cross section is a function of the top and antitop 
rapidities and 1/4 and their common transverse momentum p± . 

The triply differential cross section in ( p.5|) may be reduced to a function of a single 
variable by integrating over two of its independent degrees of freedom. We will concentrate 
upon the transverse momentum differential cross section which we obtain by integrating 
d'^a/dysdy4dp± over the rapidity range —2.5 < 2/3,2/4 < 2.5. S The resulting p± distribu- 
tion of ti pairs produced at the LHC is plotted in fig. ^. i The solid curve in the figure illus- 
trates the QCD differential cross section daQcoiPP tt)/dp±. The dot-dashed, dashed 
and dotted curves delineate the contributions from operators Og , Oj^ and O2 that are 
generated after respectively setting Cq^\A) = 0.5, C[^\A) = 0.5 and C2^\A) = 0.5 in 
Ceff with A = 2 TeV. The QCD and nonrenormalizable operator curves must be added 
together to obtain the effective field theory differential cross sections daEpriPP tt)/dp± 
that correspond to these nonzero coefficient values. 

Several points about the results displayed in fig. || should be noted. Firstly, our 
choice for the operator coefficients is simply representative. Larger coefficient values lead 
to greater differences between the QCD and EFT curves, while smaller values diminish 
the discrepancies. Our particular choice for these parameters and the scale A yields the 
combined coefficient C^^l^iA)/^^ = 0.5/(2 TeV)^ - 47r/(10 TeV)^ for the dimension-6 
operators in effective Lagrangian ( |1 . 1| ) . This value for the total coefficient is quite con- 
servative compared to that which has typically been used in previous quark substructure 



^ This rapidity range does not represent a fiducial cut but rather a reasonable integration 
interval which contains the bulk of the produced top quarks. We have checked that extending the 
range to —6 < 2/3,2/4 < 6 does not noticibly alter our final results. 

^ All of the results in fig. ^ and subsequent figures were calculated using the next-to-leading 
order parton distribution function set B of Harriman, Martin, Roberts and Stirling pT|| evaluated 
at the renormalization scale ^ = m± = \J m| + pi_'^. 
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studies. Secondly, the transverse momentum dependence of the curves in fig. ^ differenti- 
ates the dimension-6 operators from each other and the QCD terms in Ces- At low values 
of p_L, the contribution to da{PP tt)/dp±_ from the magnetic moment operator Oq^^ 

(6) 

dominates over those from the triple and double gluon field strength operators 0\ and 
0^2 . But the O^"^ curve falls off much more rapidly with increasing transverse momentum. 
So by placing a lower p± cut around 500 GeV, one can eliminate most of the chromomag- 
netic moment operator's contribution while retaining much of that from the double and 
triple field strength operators. Finally, all the curves in the figure will be shifted around 
by higher order QCD corrections. The next-to-leading 0{a^) corrections to the tree level 
QCD differential cross section may be comparable to or even larger than the 0{a1) devia- 
tions induced by the dimension-6 terms in Ces depending upon their coefficients. But the 
QCD and EFT distributions should be compared at the same order in a^. Only then can 
deviations between them be attributed to new physics beyond the Standard Model. 

It is interesting to compare the LHC differential cross sections in fig. ^ with their 
Tevatron analogues shown in fig. |[ The Tevatron curves were calculated at ^/s = 1.8 TeV 
using the same values for the nonrenormalizable operator coefficients. Not surprisingly, 
the total integrated cross-section for ti production is two orders of magnitude lower at 
the Tevatron than at the LHC. Event rate is thus more of an issue at the lower-energy 
machine. We also clearly see from the two figures that the relative importance of the 
dimension-6 terms in the effective Lagrangian depends upon collider center-of-mass energy. 
At the Tevatron, O2 dominates over Ol for equal values of their high energy scale 
coefficients. This finding is intuitively reasonable since the gluon content of colliding 
hadrons at ^/s = 1.8 TeV is less important than at ^/s = 14 TeV. 

The double gluon field strength operator affects quark scattering for all fiavors. By 
comparing its predicted impact upon the inclusive jet cross section with 1988 Tevatron 
data, we have previously set an upper bound |Cf ^(A)|/A2 < 4n/{2 TeV)^ on its coefficient 
|T|]. If Cf^ is allowed to assume this limiting value, we find that O^^ doubles the total 
integrated top quark cross section relative to the lowest order QCD prediction. This 
result is intriguing in light of the recent CDF measurement a{PP tt)cDF/o'{PP 
tt)QCD = (13.9l^:^pb)/(5.10l[; j^pb) ^ 2.7t\-l ©H- While it is premature to draw any 



conclusion from this observation, we believe it is safe to say that the bound on should 
be significantly strengthened in the future by Tevatron top quark data. 
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4. Mapping the coefficient parameter space 

In our effective field theory framework, the low energy effects of any new strong 
interaction physics are encoded into the coefficients of the nonrenormalizable terms in the 
effective Lagrangian. These coefficients define a multi-dimensional parameter space. In 
order to assess the likelihood of detecting signals of new strong sector physics at the LHC, 
we need to map this space and determine the regions where deviations from QCD could 
be measured. Within those regions, we would then like to know whether it is possible 
to isolate effects from individual operators in £efi- We will explore these issues for the 
dimension-6 operators in our basis in this section. 

To begin, we need to identify a measure of the difference between the predictions of 
QCD and the strong interaction effective theory for top quark pair production. We will 
focus upon the disparities in their LHC transverse momentum differential cross sections 
dcTQcoiPP tt)/dpj_ and daEFxiPP tt)/dp±. It is important to recall that experimen- 
tal systematic errors will inevitably render uncertain the absolute normalization for the 
observed ti distribution. In future experimental analyses, this normalization will undoubt- 
edly be fitted to QCD at low transverse momenta where any effects from new physics are 
expected to be small. We need to take this renormalization into account in our theoretical 
analysis. Therefore, we rescale by a multiplicative constant the effective theory differential 
cross section which corresponds to the distribution that will be experimentally measured. 
We choose the constant so that the renormalized daEFT / dp± cross section coincides with 
dcTQCD / dp±_ at its maximum point. 

One simple choice for a dimensionless measure of the difference between the QCD and 
EFT predictions for da{PP tt)/dp± is the ratio of their integrals: 

^ ^ fdp±{daEFT/dp±) 
J dp±{daQCD/dp±) ' 

Since the disparity between daQCD/dp±_ and daEFT / dp± increases with pj^, we perform 
the integrations in the numerator and denominator of only over the high transverse 
momentum range 500 GeV < Pjl < 1000 GeV in order to enhance the deviation of Rp^ 
from unity. The dependence of Rp^ upon the coefficient of the chromomagnetic moment 
operator Oq is then much weaker than that for gluonic operators 0\ and O2 ■ So we 
plot Rp^ as a function of C^\h) and C^\h) in fig. ^with all other operator coefficients 
set equal to zero at the scale A. 
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The origin in fig. ^ necessarily fies along the Rp^ = 1 contour, for the efi'ective field 
theory reduces to QCD at this point. It is clearly offset from the center of the concentric 
contours displayed in the figure. The offset is produced by the 0(1/A^) interference terms 
in the squared amplitude expressions ( |3.1j ) and (|3.3|) which are linear in C^^'' and €2^^ . The 
smallness of the displacement in the C| direction demonstrates that the term proportional 
to Cl ' /A^ in \A{gg dominates over the Cj /A^ interference term as we have 

previously discussed. 

While the ratio Rp^^ provides a useful global measure of the difference between the 
QCD and EFT transverse momentum tt distributions, two points in the C{ -C2 plane 
that lie along the same contour in fig. ^ may correspond to two very different da{PP — > 
ti) I dp±_ curves. It is therefore instructive to consider a second dimensionless measure 
which is sensitive to the curves' shapes. To construct such a quantity, we first discretize 
the transverse momentum interval 500 GeV <p±< 1000 GeV into = 20 bins. We then 
multiply the value for da /dp±_ in each bin by the binwidth Apj^, the integrated luminosity 
/ Cdt, the branching ratio BR for the ti pair's single lepton plus jets decay mode, and the 
6-tagging efficiency eh to convert the differential cross section into a corresponding number 
of observable ti events: 

Ni = (^) X Ap±_ X j Cdtx BR X efe. (4.2) 

We shall take the numerical values for these parameters to be Ap_i_ = 25 GeV, / Cdt = 
30 fb~^, BR = 24/81 and = 0.25. Finally, we quantify the difference between the 
effective theory distribution which will be measured for nonvanishing coefficient values 
and the theoretical QCD prediction in terms of a function. The highest p± bins are 
the most important for discriminating between the two distributions, but they contain the 
fewest events. So we adopt the Poisson function 

r A/'°'-'^ 



X 

i=l 



which is appropriate for low statistics [23]. 



We plot 'X^ jN in fig. ^ as a function of C^^''(A) and C2^''(A) over the same region 
of coefficient parameter space as in fig. ^ The innermost crescent contour in the figure 
corresponds to the expectation value x? = 1- For points lying within this contour, 
the probability that an observed transverse momentum differential cross section could be 
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attributed to a statistical fluctuation of QCD rather than to nonvanishing values for the 
Ci (A) and C2 (A) coefficients is greater than 50 %. As QCD and the strong interaction 
effective theory cannot be meaningfully told apart inside the first contour, its boundary 
establishes a limit on the sensitivity to new gluon sector physics which can be achieved at 
the LHC. The surrounding contours in fig. |^ illustrate selected "X^ jN standard deviation 
levels where a = ^J2/N = 0.316 for = 20 de grees of freedom. For example, the 
gluonic operator coefficients lying on the outermost contour yield da{PP tt) eft / dp±_ 
distributions which can be distinguished from da{PP tt)QCD / dp±_ at the 8cr level. If 
the horizontal and vertical axes were drawn to the same scale in fig. |^, this last contour 
would appear as an ellipse approximately three times more narrow in the C\ (A) direction 
than in the C2 (A) direction. The /N curves thus quantify the extent to which future 
LHC experiments will be more sensitive to the gluonic operator than to its four-quark 
competitors. 

The information contained within the contour plots of fig. ^ and fig. |^ is insufficient to 
completely determine where an observed da{PP — > tt)/dp± function lies within the C^^''- 
6*2 parameter space. But taken together, the two graphs significantly restrict the allowed 
values for these coefficients. Clearly, other differences between QCD and the effective field 
theory can be investigated along the lines which we have followed here. In particular, their 
predictions for the angular distributions of ti pairs can further constrain the allowed region 
within the coefficient parameter space. We briefiy touch on this topic in the next section. 



5. Top-antitop angular distributions 

Our study of the prospects for detecting new strong interaction physics at the LHC has 
so far utilized only the transverse momentum information incorporated within d'^a{PP —>■ 
ti) /dy^dy^dp^. The triply differential cross section contains, however, complementary 
angular distribution information. We consider its implications for discriminating between 
the nonrenormalizable operators within the effective Lagrangian in this section. 

A number of angular distributions for ti pairs as well as their decay products can 
be generated by integrating d^a / dy3dy4^dp±_ over various ranges of 1/3, 1/4 and p_L. One 
differential cross section of particular interest is da{PP ti)/d cos 9* where 9* denotes 
the angle between the direction of the boost and that of the top quark in the parton center- 
of-mass frame. We plot this cross section in fig. |^ for pure QCD and QCD plus some of the 
nonrenormalizable operators in our basis set. In order to enhance the operators' signal over 
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the QCD background, we have imposed the transverse momentum cut p± > 500 GeV. We 
have also required the lab frame angle between the tort and the beamline to exceed 25.4° 
to approximate the acceptance of planned LHC detectors. This last restriction ensures that 
the pseudorapidities of the decay products from high momentum tops will predominantly 
lie within the interval —2.5 < r] < 2.5. 

The solid curve in fig. ^ represents the QCD differential cross section da{PP — > 
tt)QCD / d cos 6* . The dotted, dashed and dot-dashed curves in the figure illustrate 
da{PP tt)„pT/dcose* for cf^(A) = 0.5, Cf^(A) = 0.5, and Cf^(A) = -0.5 with 
A = 2 TeV. We again find that the effect of the triple gluon field strength operator 

(6) 

Oi is essentially independent of the sign of its coefficient. The differential cross section 
corresponding to c[^\A) = —0.5 thus closely traces that for C'^^)(A) = 0.5 displayed in 
the figure. The dimension-8 gluon operator O3 induces deviations from pure QCD which 
are clearly visible in da /d cos 9*. This result is both surprising and interesting since we 
previously found that the effect of O3 upon the tt transverse momentum distribution 
was negligible. Indeed, we did not include a curve corresponding to the dimension-8 gluon 
operator in the transverse momentum plots of fig. |^ and fig. ^ since it would have been sup- 
pressed relative to its dimension-6 counterpart by more than an order of magnitude. But 
we now see that experiments at the LHC can be sensitive to this next-to-next-to-leading 
order operator if they probe its impact on the ti angular distribution. 

As for the transverse momentum distributions, it is again useful to identify di- 
mensionless measures of the differences between the QCD and EFT predictions for 
da{PP tt)/d cos 9*. One simple choice for such a measure is the ratio of their inte- 
grals 

/ d cos 9* {da E FT /d cos 9*) 
" J d cos 9* {daQco/d cos 9*) ^ ' 



which is the analog of Rp^ in eqn. Another is the ratio i?rms of their root- 

mean-squared values for cos^*. We have calculated these ratios for the angular differ- 
ential cross sections shown in fig. ^ and for similar cross sections involving the mag- 
netic moment and four-quark operators. For the curves corresponding to Cq (A), 
Cf^(A), Cf^(A) and Cf^(A) equal to 0.5, we find i?ang = (1.03,1.23,1.46,0.82) and 
-Rrms = (0.999,0.978,0.991,0.871). For analogous curves with the A scale coefficients set 
equal to -0.5, we find R^^g = (0.969, 1.23, 0.735, 1.18) and i^^ms = (1-00, 0.978, 1.03, 1.08). 



7 (8) 

The coefficient Cg was not evolved using the renormalization group but was instead simply 



^3 

fixed at its A scale value. 
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The most striking conclusion which we draw from these results is that the dimension-8 
gluon operator alters the shape of the tt angular distribution much more than all the other 
dimension-6 operators in >Ceff for comparable values of their coefficients. The magnetic 
moment operator's angular distribution is indistinguishable from that of pure QCD, while 
the distributions of the and {DGY operators differ significantly from that of QCD in 
-Rang but not in -Rrms- O^^^ thus possesses a distinctive signature: a QCD- like transverse 
momentum distribution but a quite nonstandard angular distribution with -Rang and -Rrms 
both deviating from unity in the same direction. 

A more precise determination of the power of angular measurements to delimit allowed 
regions of coefficient parameter space will require detailed simulations including top quark 
decays and detector resolution. We leave such a study to future work. 

6. Conclusions 

In this article, we have investigated the impact of the triple gluon field strength 
operator upon top quark pair production. The operator represents the only genuinely 
gluonic CP even term which can arise at dimension-6 within an effective strong interaction 
Lagrangian. Although it has proven surprisingly difficult to study in the past, the prospects 
for either detecting or significantly constraining this operator appear quite promising at 
the LHC where its effects will be significantly enhanced by the large gluon content at 
small X of the colliding hadrons. We have found that the sensitivity of the tt transverse 
momentum distribution to the triple field strength operator is greater than that to all other 
dimension-6 pure quark and mixed quark-gluon terms in >Cefr for comparable values of their 
coefficients. Moreover, only one higher order gluonic operator can contribute to gg — > tt 
at dimension-8, and its effect on the distribution is more than an order of magnitude 
smaller than that of G^. We have also seen that angular distribution information can help 
to differentiate effects from the various operators which may reside within the effective 
Lagrangian. Top-antitop production therefore promises to provide an important means 
for probing new strong interaction physics beyond the Standard Model. 
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Figure Captions 



Fig. 1. Diagrams witli a lieavy colored fermion or boson running around the loop which 
match onto gluonic operators in the effective strong interaction Lagrangian. 

Fig. 2. (a) Lowest order QCD graphs which contribute to gg ti scattering, (b) 
0(1/A^) gg —>■ ti graphs with single insertions of the chromomagnetic moment 
operator oj^\ (c) 0(1/A^) gg ^ ti graph with single insertion of the triple gluon 
field strength operator o\^\ (d) 0(1/A^) gg ti graphs with two insertions of 
either oj^^ or O^^^ and 0[^\ 

Fig. 3. Lowest order QCD, chromomagnetic moment and four-quark operator graphs 
which contribute to qq ti scattering. 

Fig. 4. Transverse momentum differential cross section da{PP ti)/dp± calculated for 
an LHC center-of-mass energy -^i = 14 TeV. The solid curve illustrates the pure 
QCD cross section. The dot-dashed, dashed and dotted curves represent the 
additional nonrenormalizable operator contributions obtained after respectively 
setting ci^^(A), cf ^(A) and ^(A) equal to 0.5 with A = 2 TeV. 

Fig. 5. Transverse momentum differential cross section da{PP ti)/dp± calculated for 
a Tevatron center-of-mass energy ^/s = 1.8 TeV. The curves are labeled the same 
as those in fig. ^. The dashed O^^ curve must be subtracted from rather than 
added to the solid QCD curve to obtain the EFT cross section that corresponds 
to cf^(A) = 0.5. 

Fig. 6. Ratio Rp^ of the EFT and QCD ti transverse momentum distributions integrated 
over 500 GeV < p± < 1000 GeV plotted as a function of gluonic operator coef- 
ficients c[^\A) and ^(A) with A = 2 TeV and = 14 TeV. The values for 
Rp^ are displayed alongside the contours. 

Fig. 7. x^/A^ for = 20 transverse momentum bins plotted as a function of gluonic 
operator coefficients cf ^ (A) and ^ (A) with A = 2 TeV and = 14 TeV. 
The innermost crescent contour corresponds to x^/A^ = 1. The surrounding 
contours represent /N = 1 + 2ncr where a = = 0.316 and n = 1,2, 3, 4. 

Fig. 8. Angular differential cross section da{PP ti)/d cos 9* calculated for y/s = 14 
TeV. The solid curve illustrates the pure QCD cross section. The dotted curve 
corresponds to QCD plus the dimension-6 operator O^^ with c[^\A) = 0.5 
and A = 2 TeV. The dashed and dot-dashed curves represent QCD plus the 
dimension-8 operator O3 with C3 (A) set equal to 0.5 and -0.5 respectively. 
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